Absorption/Expulsion of Oligomers and Linear Macromolecules in a Polymer Brush 
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The absorption of free linear chains in a polymer brush was studied with respect to chain size 
L and compatibility \ with the brush by means of Monte Carlo (MC) simulations and Density 
Functional Theory (DFT) / Self-Consistent Field Theory (SCFT) at both moderate, a g = 0.25, and 
high, o~ g — 1.00, grafting densities using a bead-spring model. Different concentrations of the free 
chains 0.0625 < (f> < 0.375 are examined. 

Contrary to the case of \ = when all species are almost completely ejected by the polymer brush 
irrespective of their length L, for \ < we find that the degree of absorption (absorbed amount) 
r(L) undergoes a sharp crossover from weak to strong (« 100%) absorption, discriminating between 
oligomers, 1 < L < 8, and longer chains. For a moderately dense brush, a g = 0.25, the longer species, 
L > 8, populate predominantly the deep inner part of the brush whereas in a dense brush a g — 1.00 
they penetrate into the "fluffy" tail of the dense brush only. Gyration radius R g and end-to-end 
distance R e of absorbed chains thereby scale with length L as free polymers in the bulk. Using both 
MC and DFT/SCFT methods for brushes of different chain length 32 < TV < 256, we demonstrate 
the existence of unique critical value of compatibility x ~ X° < 0. For X c (4'o) the energy of free 
chains attains the same value, irrespective of length L whereas the entropy of free chain displays a 
pronounced minimum. At \ c all density profiles of absorbing chains with different L intersect at 
£J5 ' the same distance from the grafting plane. 

The penetration/expulsion kinetics of free chains into the polymer brush after an instantaneous 
change in their compatibility \ displays a rather rich behavior. We find three distinct regimes of 
penetration kinetics of free chains regarding the length L: I (1 < L < 8), II (8 < L < TV), and III 
(L > TV), in which the time of absorption r grows with L at a different rate. During the initial stages 
of penetration into the brush one observes a power-law increase of V oc t a with power a a — In <f) 
whereby penetration of the free chains into the brush gets slower as their concentration rises. 
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I. INTRODUCTION 

> 

(T) ', Densely-grafted chains on nonadsorbing substrate surfaces form the so-called "polymer brush" [Il4l6| . These systems 
find various important applications [i"6j . e.g. as lubricants [HJ], for colloid stabilization [l7|, for tuning of adhesion and 
wetting properties [l6l 1 1 8| , for improving the biocompatibility of drugs [l9j , as protective coatings preventing protein 
adsorption ( "nonfouling" surfaces) in a biological milieu [20|, microfluidic chips for biomolecule separation [21], etc. 

The theoretical description of the conformations of macromolecules in these polymer brushes and their dynamics 
has been an active topic of research hitherto (e. g., l22|-[28l| ; for early reviews see |9Hl5j. Also the interaction of the 
brushes with either the solvent molecules (e.g. [2^1 j25i Tor globular proteins [24| and/or other nanoparticles (e.g., 
J> ! [HI HE I29rl33j ) has found much recent attention. However, in many situations of interest there will also occur free 

k>( ' polymer chains in the solution, interacting with the polymers of the brush. This interaction has received relatively 
less attention, apart from the case where a polymer brush interacts with a dense polymer melt @, [33443} . The latter 
case is particularly interesting because there is very little interpenetration of the grafted chains in the melt and the 
free chains in the brush even if their chemical nature is identical ( "wetting autophobicity" (i3 - |47j ) . 



In contrast, scaling theory [2j, self-consistent field [48[ and simulation [49| have predicted partial penetration of 
free chains into moderately dense brushes of identical chemical nature in semi-dilute solutions when the monomer 
volume fraction in solution approaches that of the brush, and this behavior has been confirmed experimentally [4a |. 
Of course, when the polymer solution is very dilute, the brush provides a free energy barrier for penetration of free 
chains into it and this limits the grafting density that can be achieved when one prepares a brush by grafting chains 
from solution [la ] (see also some attempts to model this process by simulations [50L l5l|). Similarly, since typically 
the energy won by the chain end when it gets grafted is of the order of lOkrjT only jldl . [l2l [ill ], there is a nonzero 
probability that brush chains get released from the grafting substrate surface and are subsequently expelled from the 
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However, most cases studied so far refer to the situation that (apart from chain end effects) the chains in the bulk 
and those in the solution are identical. It is interesting, therefore, to consider the more general situation when the 
grafted chains and those in the bulk differ in their chemical nature. Then the problem of compatibility (traditionally 
modelled by introducing a Flory-Huggins ^-parameter |53|) between the two types of chains arises. Then, there is 
also no reason to assume that the length N of the grafted chains, and the length L of the free chains are equal. 
Such situations (in particular, when the grafted and the free chains attract each other, \ < 0) are of great interest 
for modern applications such as protein adsorption, " antifouling " surfaces [20], etc. However, to the best of our 
knowledge, no systematic study of the effects of the various parameters (N, L, \ and monomer concentration of the 
free chains </> Q ) on the amount of absorption and the penetration kinetics has been reported so far. The present paper 
presents simulation and Density Functional Theory (DFT) results in an effort to fill this gap. In Section|H]we describe 
the model and comment on some simulation aspects; Section [IIII summarizes our theoretical approach which includes 
both static and dynamic versions (DDFT) of DFT as well as Self-Consistent Field Theory (SCFT) calculations. The 
numerical results are described in Section llVl while Section \V\ contains a summary and discussion. 



FIG. 1: Snapshots of a polymer brush with chain length TV = 32 at grafting density a g — 0.25 and free chains of length L — 32 
at equilibrium: (left) absorption at e po = 2.00, and (right) expulsion at e po = 0.01. 



We have used a coarse grained off-lattice bead spring model[55l. |56j to describe the polymer chains in our system. 
As far as for many applications in a biological context rather short grafted chains are used (57j . we restrict ourselves 
to length N = 32 of the grafted chains. The polymer brush consists of linear chains of length N grafted at one end to 
a flat structureless surface. The effective bonded interaction is described by the FENE (finitely extensible nonlinear 
elastic) potential, 



with K = 20, Imax — 1, lo = 0.7, l m in = 0.4. Thus the equilibrium bond length between nearest neighbor monomers 
is lo — 0.7. Here and in what follows we use the maximal extension of the bonds, l m ax = 1, as our unit length while 
the potential strength is measured in units of thermal energy ksT where fc^ is the Boltzmann constant. 
The nonbonded interactions between brush and free chain segments are described by the Morse potential, 




II. MODEL AND SIMULATION ASPECTS 
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with a = 24, r m i n = 0.8, andtM /ksT standing for the strength of brush-brush, e pp , polymer-polymer, e oa , and 
brush-polymer, e po interactions. In our present study we take typically e pp = 0.2, e oa — 0.1 (that is, in the good 
solvent regime with only excluded volume interactions). For e po = 2.00, \ = —1.85 the free chains are absorbed in the 
brush whereas for e po — 0.01, \ — —0.005 the polymer brush ejects them into the bulk. Note that we define here the 
compatibility parameter \ simply as y = 0.5(e pp + e 00 ) — e po , and do not include the coordination number (which is 
done when one uses lattice models [531 ] . 

The size of the container is 16 x 16 x 32. The polymer chains are tethered to grafting sites which constitute a 
triangular periodic lattice on the substrate whereby the closest distance between grafting sites is l m ax — 1- Thus 
the largest grafting density a g — 1.0 involves 8192 brush segments, if the polymer chains are anchored at distance 
lmax = 1) and o g = 0.25, i.e., 2048 brush segments, if the 'lattice constant', i.e., the distance between adjacent head 
monomers on the surface is equal to 2l max . Note that a = 1.0 corresponds to a simulation where the monomer density 
in the brush near the wall is close to the density of a polymer melt while a = 0.25 would correspond to a rather 
concentrated polymer solution. 

For the chain model, em/hsT = e pp — 0.2 corresponds to good solvent conditions since the Theta-point for a 
(dilute) solution of polymers described by the model, Eqs. [l][2]has been estimated [55l] as &b©Am = 0.62. In all our 
simulations we use brushes formed by polymer chains consisting of N = 32 effective monomers whereas the number 
of free chains N a of length L (where L spans the interval 1 < L < 64) is taken such that the total number of free 
chain segments remains constant and is equal to 512. For a certain length L = 64, however, we also change the 
concentration of free chains in the container by varying their number N a in the interval 8 < N D < 48. Thus, the 
volume fraction c G of 64-free chains is varied between 0.0625 < c < 0.5. Note that, as usual, solvent molecules are 
not explicitly included (56l . l58l . [59[ but work which includes solvent explicitly [23| would yield very similar results. 

For a dense brush with polymer chains of lengths N = 32 statistical averages were derived from typically 10 7 Monte 
Carlo Steps (MCS) per monomer. The Monte Carlo algorithm consists of attempted moves whereby a monomer is 
chosen at random and one attempts to displace it to a new randomly chosen position —0.5 < Ace, Ay, Az < 0.5 
regarding the old position. We use periodic boundary conditions in the x — y directions and impenetrable hard walls 
in the z direction. Two typical configurations of the polymer brush with free chains of length L = 32, are shown in 
Fig. Q]for the case of good, e po = 2.00, and poor, e po = 0.01, compatibility with the polymer brush. 



III. THEORY 



We employ classical DFT to compute density profiles of free and grafted polymer chains. Theory has been discussed 
in detail in previous publications, so here we briefly summarize its most important aspects. The starting point of 
the DFT treatment is the expression for the grand free energy, J7, as a functional of the density profiles of free and 
grafted chains, </> (R ) and </> p (R p ), respectively (R p / = {r,}, where are the positions of the chain segments). The 
functional f2 is related to the Helmholtz free energy functional, F, via a Legendre transform: [60l l6~i| 

f>[0 o (R o ),0 p (R p )] = J F[0 o (R o ),0 p (R p )]+ I dRMR a )V a (R a ), (3) 

a—o,p 

where V a (R a ) is the external field, which in the present case is due to the hard-sphere like interaction of the polymer 
segments with the hard wall, Vp(R p ) = Y]^-i v p (zi), where v p (zi) = oo for z < and v p (zi) = otherwise, with 
analogous expression holding for V (R ). Additionally, the innermost (i = 1) bead of each grafted chain is tethered 
to the wall via a grafting potential exp[— /3w p (zi)/fcsT] = 5(zi), where j3 — l/fc^T. Note that the chemical potential 
of both free and grafted chains is absent from the second term of Eq. §3§ because the DFT calculations are performed 
at a fixed number of both free and grafted segments in order to mimic the MC simulations: / "* ax dz4> p (z) = a g N 
and / Zmax dz<j) (z) = N D L/A. in the above, z max — 32 is the box length and A = 256 is the wall area. 

The Helmholtz free energy functional is separated into ideal and excess parts. [60l l6l| with the former given by: 

/3F ld [0 o (R o ),0 p (R p )] = { / dRa0a(Ra)[ln0 Q (Ra) - 1] + P / dR a <p a (R a ) V b (K a ) I . (4) 

where the bonding energy Vb for the grafted chains is taken as follows: 

r mi in w i TT 1 5 (l r i ^ r i+i| — ^) fc -s 
cxp[-f3V b {Rg)} = [[ — j2 , (5) 

with a similar expression for free chains, with N replaced by L. This bonding potential constrains adjacent segments 
to a fixed separation 6/. 
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The excess part of the Helmholtz free energy is written as a sum of repulsive (hard chain) and attractive terms, 
with the former computed in the weighted density approximation and the latter obtained within mean-field approach, 
using Eqs. (12)-(17) from Ref. [Ill; for the sake of brevity we do not reproduce these equations here. 

The minimization of the grand free energy functional with respect to P (R P ) yields the equilibrium density distribu- 
tion for the grafted chains which can be integrated over grafting and bonding delta-functions to obtain the following 
result for the density profile of the ith segment of the grafted chains: [13] 

cf> pi (z) = C i I p (z)Ir(z)lf(z), (6) 

where 

I p (z) = exp[-p(v p (z) + \ p (z))l (7) 

with 

b(p p {z) 

The two propagators in Eq. ©, if and I~ move from the free (i = N) and the tethered (i = 1) ends of the chain, 
respectively. They are computed via recursive relations given by Eqs. (23)-(25) of Ref. [Hj]. 

The normalization constant Cj in Eq. is chosen to ensure that the ith segment density profile is normalized to 
cr g . The total segment density profile for the grafted chains is given by: 

N 



(z) =J2<f> P i(z)- (9) 



The equilibrium density profile for the segments of the free chains can be obtained in a similar way, by minimizing the 
grand free energy functional with respect to O (R O ) and integrating over bond- length constraining delta functions. 

The DFT equations described above are solved simultaneously to obtain the segment density profiles for free 
and grafted chains. The equations are solved iteratively using Picard algorithm, [25| with the step size along the 
z coordinate taken to be 0.0325. The above procedure yields equilibrium segment density profiles for a given set 
of interaction potentials. In addition to the equilibrium structural properties, we have also studied the kinetics of 
the adsorption of free chains into the brush, following a switch of the interaction potential between free and grafted 
segments from repulsive to attractive. To this end, we have employed the DDFT method, which is a dynamical 
generalization of the DFT approach. (gH |(54| MC simulations have indicated that the segment density profiles of 
the grafted chains are essentially independent of the strength of the attraction between free and grafted segments. 
Accordingly, in our DDFT calculations we take <f> p {z) to be time independent and focus on the time dependence of 
the free chain density, (j> (z,t). 

The time evolution of the segment density profile of free chains is given by the following equation: [64j 

where ^,{z,t) is the non-equilibrium local chemical potential, and dimensionless time r is defined according to r = 
ksTM/l^^t, where M is the mobility coefficient. 

Initial density profile <j> {z,t = 0) corresponds to the equilibrium distribution of free chains at a repulsive brush, 
i.e. e po = 0.01. At t — 0, the brush-free polymer attraction is instantaneously "switched on", i.e. we set e po = 2. 
The time-dependent polymer density profile is then propagated according to the Eq. (jlO[) . with the time-dependent 
chemical potential given by: 

L 

0li(z, t) = In 4> (z, t)-\nJ2 Cil {z, t)/r (z, t)I+(z, t), (11) 

i=l 

where I {z, t) is obtained by substituting the time-dependent density 4> (z, t) into the expression for I (z) (and likewise 
for the propagators I~{z 7 1) and lf(z, t). We solve Eq. (ITU1) using Crank-Nicholson scheme. [63ll64{ Note that Eq. (fT0|) 
has the form of a continuity equation with the flux (current density) given by j(z,t) — —(f> (z,t)^-/3p l (z,t). The fact 
that the DDFT method propagates 4> {z,t) via a continuity- type equation guarantees the conservation of the total 
number of segments in the system, which is consistent with the simulation set-up. 

In order to compare the results of the DDFT approach with kinetic MC data, we set the mobility coefficient M 
equal to unity and adjust the conversion factor between DDFT dimensionless time and kinetic MC number of steps 
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for one particular set of parameters L and <fi . Comparisons for all other values of L and <f> are performed using the 
same conversion factor, while assuming M to be inversely proportional to both L and <p . 

With the goal of shedding further light on the thermodynamic aspects of the adsorption process, we have also 
performed self-consistent held theory (SCFT) calculations of the structural properties as a function of the interaction 
strength between the segments of the brush and the free chains (in SCFT approach this interaction is characterized 
by the parameter \ which is calculated in the standard fashion from the corresponding potential well-depths: x — 
0.5 * (e pp + e 00 ) — e po ). The main motivation behind carrying out SCFT calculations is the fact that this approach 
provides a more straightforward way to decompose the free energy into entropic and energetic components, thereby 
providing a complementary (to DFT) view of the adsorption process. 

The basic equations of the SCFT method are well known, [42|, EH and will not be reproduced here for the sake of 
brevity. Once again, the density profiles for free and grafted chains are written in terms of the propagators, the only 
major difference from the DFT approach being that instead of the equation of state, one employs the incompressibility 
constraint to set up the equations for the density profiles, which are once again solved iteratively using Picard's method. 
For example, the equation for the density profile of the grafted chain segments takes the form: 



Mz) = c P }^ 1 w:: , (12) 

where the normalization constant C p is obtained from the grafting density er ff , and G p {z) = exp(—/3u p (z)), with 
j3u p (z) = u'{z) + x po ^> (z). The hard core potential u' is independent of the segment type and serves as a Lagrange 
multiplier enforcing the incompressibility condition, meaning that the lattice space is completely filled and no segment 
overlap occurs. The density profile of the free chain segments is obtained in a similar way. 

Once the profiles are calculated, one can easily obtain excess entropy and energy of the free chains (relative to pure 
unmixed components) as follows [65[: 

S-S* = -k B f maX dz^(2)(^+lnG (z)l, (13) 



U-U* = -k B T I dzMz)x P oM z ), ( 14 ) 
where <b b n is the bulk volume fraction of free chains. 



IV. RESULTS 



A. Equilibrium properties 

In Fig. [2] we show the density profiles of the free chains, <f> (z), of length L for an attractive, e po = 2.00, and a 
neutral, e po — 0.04, brush along with the monomer density profile of the brush itself, 4> p (z). Our MC simulation results 
indicate that at fixed segment concentration, Cj, the brush profile, <f> p (z), is virtually insensitive to L, whereupon we 
keep only one such profile in the graphs. The most striking feature which may be concluded from Fig. [2] is, somewhat 
counter- intuitively, the strong increase of absorption with growing length of the absorbed free chains L. Evidently, 
both at moderate, a g — 0.25, and high, a g = 1.00, grafting density, the longer polymers are entirely placed inside the 
polymer brush whereas the much more mobile short species L < 8 remain uniformly distributed in the bulk above 
the brush end. Since some of the absorbing chains with larger L get stuck inside the brush, their density profiles 
could not smoothen sufficiently for the time of the simulation run. Therefore, we observe rather large statistical 
fluctuations in (f> (z). For repulsive brushes all species are largely expelled from the brush whereby the situation is 
reversed as far as the free chain length L is concerned. In the very dense brush a g = 1.0, the brush profile displays the 
typical oscillations near the grafting surface suggesting some layering immediately in the vicinity of the grafting wall 
- Fig. [5J:,d. In all graphs one observes pronounced depletion effects at the upper container wall, opposing the brush. 
However, the inhomogeneity of 4> {z) near the wall at z — 32 has no effect on (f> (z) in the region of the polymer 
brush, the flat part of 4>o{z) in between the brush and the confining wall at z — 32 is broad enough to eliminate any 
finite-size effects associated with the finite linear dimension of the simulation box in z~ direction. 

One should note also the good agreement between simulation and DFT results. In fact, the thin lines, indicating 
the latter, may hardly be distinguished from the Monte Carlo data (thick lines) in Fig. [5^,,b. The only significant 
discrepancy between theory and simulation is observed in the brush profile in the vicinity of the grafting wall, where 
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FIG. 2: Density profiles of the polymer brush, 4> p (z), (shaded area) and of free chains, (f>o{z), (thick lines) of length L, (given 
as parameter) at two grafting densities: a g = 0.25 (upper row), and a g = 1.00 (lower row), (a) and (c) illustrate good 
compatibility between brush and free chains, e po = 2.0 while (b) and (d) demonstrate a case of bad compatibility, e po = 0.04. 
Thin solid lines in (a) and (b) denote results from the DFT calculation. The densities in (a) are normalized so as to reproduce 
the correct ratio of brush to free chains concentrations <j)p an d 0o (the absolute particle concentration a is indicated in the 
alternative y— axis. For the sake of better visibility, in (b), (c), and (d) the density of all species is normalized to unit area. 



DFT approach overestimates the oscillations. This discrepancy is likely due to the fact that in the DFT method the 
bond lengths are constraint via delta-functions to a constant value of bi = 0.75, while in the simulations the bonds are 
allowed to vibrate under FENE potential, Eq. ([1]). For L = 64, Fig. [3] shows a qualitatively similar behavior of the 
density profiles for the cases of gradually increasing free chain concentration (indicated by the number of free chains 
Af as parameter). Expectedly, for N > 24 (which corresponds to monomer concentration a — 0.1875) and a g — 0.25, 
the free chains are present in the bulk over the brush as the brush interior is then entirely filled. However, when the 
brush - free chain attraction increases to e po = 3.00, the MC data (not shown here) indicate complete absorption of 
the free chains into the brush with virtually no free chains in the bulk above the polymer brush even at the highest 
concentration of Cj = 0.375. 

With increasing grafting density and/or free chain concentration, the agreement between DFT and MC deteriorates 
somewhat, with the theory underestimating the degree of penetration of free chains into the brush (see discussion of 
Fig. [3] below), which is likely due to the simple Tarazona's weighting function employed in our DFT calculations. It 
is well known that at higher densities it would be more appropriate to use weighting functions from the Fundamental 
Measure Theory. [6^] Indeed, precisely such approach has been recently used to study adsorption and retention of 



spherical particles in polymer brushes 54 1 



Next, we present MC and DFT results for the absorbed amount of free chains as a function of degree of poly- 
merization and concentration. The absolute absorbed amount is defined as the number of polymer segments located 
"inside the brush" , namely, in the region z < z cut , where the cutoff distance z cut is defined in such a way that 99% 
of the brush segments are located in the region z < z cut . The relative absorbed amount T is defined as the ratio 
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FIG. 3: Density profiles of the polymer brush, 4> p {z), (shaded area) and of free chains, <j)o( z ), (thick lines) of length L — 64 for 
different free chain concentration (number of free chains N c ) (given as parameter) at two grafting densities: a g = 0.25 (upper 
row), and a g = 1.00 (lower row). Thin lines in (a) denote DFT results. The inset in (b) shows the change in the density profile 
of brush monomers for two values of free chain concentration: a = 0.00, and 0.375. In (a), (c) and (d) <f> p stays practically 
constant. 



of the absolute absorbed amount to the total number of free chain segments. In Fig. one may observe the steep 
increase in T with growing polymer length L both for brushes with a g — 0.25 and a g = 1.00 when e po = 2.0. Indeed, 
as indicated also in Fig. [3J as soon as L > 8, the adsorbed amount saturates at nearly 90%. A much more gradual 
growth of T is found for the critical attraction e po = 1.70 (see below). In Fig. 0^ one sees again that DFT results 
for the absorbed amount of polymers as a function of the absorbate polymerization index (shown here for the case of 
lower grafting density) are in good agreement with MC data, with the exception of the intermediate-length chains, 
where DFT overestimates the adsorbed amount somewhat. 

Especially interesting is the observation, Fig. ^p, that the conformations of the absorbed chains inside the brush 
practically do not change with respect to those of the free chains in the bulk - the scaling behavior of the parallel and 
perpendicular components of the end-to-end (squared) distance R\ and radius of gyration, R 2 g) is demonstrated in 
logarithmic coordinates by straight lines whereby the value of the Flory exponent v w 0.64. Due to the short lengths 
of the free chains used here this value is slightly larger than what is expected for very long chains (namely v = 0.59). 
Only the absorbed chains that are longer than the polymers of the brush, L = 64 > N = 32, indicate deviations from 
the scaling law of single polymers with excluded-volume interactions: the parallel component R gX y slightly exceeds, 
and the perpendicular component, R gz , falls below the straight line suggesting that the original shape of the L = 64 
coil flattens parallel to the grafting plane. 

Fig. [5] displays the dependence of absorbed amount of polymers on the concentration for the highest polymerization 
index studied, L = 64. One sees that for both grafting densities the total number of absorbed monomers increases 
with concentration, while the relative absorbed amount decreases. DFT results (again presented for the case of lower 
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FIG. 4: (a) Variation of the absorbed amount F with polymerization index L of the free chains for two grafting densities. Empty 
symbols denote DFT results. The case e po = 1.70 refers to the critical degree of brush-polymer compatibility (cf. Section TlV B|) . 
(b) Mean squared radius of gyration, R gxy , and end-to-end distance, R^. xy , parallel and perpendicular, R gz , R\ z to the grafting 
plane against length of the free chains L at grafting density a g = 0.25. Dashed lines denote the observed slope ~ 1.28 ± 0.03. 
Only for the longest free chains with L = 64 a marked deviation from the standard scaling behavior may be detected. 
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FIG. 5: Variation of the absorbed amount T with free chain concentration c; for L = 64 and two grafting densities: a g = 0.25 
(circles), and a g = 1.00 (squares). Full symbols correspond to polymer absorption with e po = 2.0 and empty symbols denote 
expulsion e po = 0.01. (a) Absorbed fraction vs. a, (b) Total number of absorbed monomers against d. 



grafting density) fall below MC data at higher concentrations, illustrating the aforementioned observation that DFT 
underestimates the degree of penetration of free chains into the brush at higher concentrations. 



B. The critical compatibility \ c 



As a remarkable feature of polymer absorption in a brush we find the existence of a critical degree of compatibility 
X° between the grafted and free chains. Fig.[SJi displays brush and free chain density profiles for various polymer chain 
lengths at the critical value of the brush-polymer attraction strength [e^ = 1.70 for MC simulations and Xpo — —1-40 
for SCFT calculations). While simulation and theoretical results differ quantitatively, there is a striking qualitative 
similarity in that the density profiles, irrespective of the length L of the free chains, all intersect in two single points 
(inside and outside the brush). The DFT approach produces exactly the same behavior albeit for a smaller fl = 1.0 
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FIG. 6: (a) MC Density profiles 4>o{z) at the "critical" strength of attraction e c vo = 1.70 for different lengths L. (b) SCFT 
results for the variation of energy (solid lines) and entropy AS (symbols) of free chains of length L with changing attraction 
X P o to the polymer brush. Arrow indicates the intersection point of energy, Xpo = —1-40, which coincides with the position of 
the minima in AS. All energy values are multiplied by 10 for better visibility. In the inset the entropy AS for chains with 
L = 64 in the brush (full squares) and in the bulk (empty squares) is displayed against e po . 



(not shown here). Fig. [BJd shows SCFT results for the excess entropy and for the internal energy per monomer (given 
by Eqs. (IT31 and (|1~4"]) . respectively) as a function of Xpo- One notes immediately that all the energy curves intersect 
in a single point, corresponding to Xpo = Xpoi while all the entropy curves pass through a minimum at this point. 
Furthermore, the entropic curves corresponding to the polymer segments located "inside" and "outside" the brush 
(as defined earlier) intersect at the same value of Xpo as shown in the inset of Fig. |6Jd. 

While at Xpo there exists thus a distance z from the grafting plane where the local concentration of polymer 
solutions is independent of polymer length L, provided <p po is kept constant for all L, the value of x%o itself is expected 
to depend on the concentration and/or the size of the grafted chains N. We performed SCFT calculations to see how 
the "critical" value of Xpo changes with a and N within a broad range: < a < 0.375 and 32 < N < 256. We 
find that it increases as Xpo = 1-306 + 1.326c, + 2.393cf with increasing free chain concentration a, and decreases as 
Xp = 1.874iV -0 0858 with increasing length N of the grafted chains (in the latter case, the grafting density is adjusted 

1/3 

such that the typical scaling variable for grafted polymers Na g ' is kept constant). 



C. Adsorption/Desorption Kinetics 



Here we present our simulation and theoretical results for the kinetics of polymer adsorption/desorption into, or 
out of the brush. Fig. [7^, shows the variation of the absorbed relative amount, T, with elapsed time t following an 
instantaneous switch of the interaction between brush and free chains. As expected, the expulsion of the adsorbate 
from the brush after an instantaneous switching off of brush - polymer attraction proceeds much faster than the 
absorption kinetics. The latter, as is visible from the inset to Fig. [7^,, proceeds through an initial steep increase toward 
a saturation plateau of T whereby the small species absorb faster than those with larger L. From the intersection of 
the tangent to the initial steep growth of T and the saturation value one may determine the characteristic time of 
absorption r as function of L - Fig. [JJd. The results are presented for all values of free chain lengths, and one sees 
that DDFT results are again in good agreement with kinetic MC data. This also holds in Fig. [7t> where indeed the 
theory is in good agreement with simulations for N = 32. For the case of longer grafted chains (iV=64, <r 9 =0.2), 
no simulations were performed and only theoretical predictions are shown. Nonwithstanding, for both values of N, 
one clearly sees three regimes in the dependence of r on L. In the first regime, the absorption time grows fast and 
essentially linearly with L (up to L = 8 for the shorter brush and L = 16 for the longer one). By analyzing the data 
presented in Fig. |4Jd, one sees that this initial linear regime corresponds to the situation when R g of absorbed chains 
is less than or equal to the average distance between the grafting points. As L (and, consequently, R g ) is increased 
beyond the aforementioned values, one enters the second regime where the growth of r, while still nearly linear, is 
markedly slower. We interpret this slowing down as a halmark of an increased friction of the penetrating coils when 
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FIG. 7: (a) Variation of the absorbed amount T with elapsed time t after an instantaneous change of the interaction between 
brush and free chains. Here a g = 0.25 and the averaging was performed over 50 cycles. The inset shows the filling kinetics for 
different size L of free chains. The total number of free chain monomers 512 was kept constant, (b) Absorption time r against 
polymer length L for a g = 0.25 displays three distinct regimes I (1 < L < 8), II (8 < L < N, and III (L > N) (shaded areas). 
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FIG. 8: (a) MC data for the absorbed amount F against elapsed time t after the onset of absorption for different concentration 
of free chains with L — 64. The log- log plot shows that V grows by power law F(t) oc t a . The measured slopes a are plotted in 
the inset against the number of free chains N . One finds a oc — | ln<^ . (b) The same as in (a) but at the "critical" attraction 
£ po — 1.70 and fixed Cj = 0.0675 where T(t) oc t 13 . The exponent /3 oc L 1 / 3 (inset). 



their radius of gyration exceeds the size of the cavities in the polymer brush. This regime extends up to the point 
where the lengths of free and grafted chains become equal. Beyond this point, for L > N, the third regime is sets in, 
where the absorption time is essentially independent of the free chain length. One might see therein an indication of 
a change in the mechanism of free chain penetration into the brush with thickness H < R g whereby additionally the 
coil flattens inside the grafted layer due to gain in absorption energy. 

Fig. [SK displays simulation and theoretical results for the absorption kinetics for N = 32, L = 64, and several 
values of the concentration Ci = 647V D /8192. Both MC and DDFT data show that at early and intermediate times the 
time dependence of the absorbed amount follows a power law T(t) oc t a . The corresponding effective exponent a is 
decreasing as the concentration increases (see inset), although the value of T at the beginning of the intermediate time 
regime is larger for larger values of N . This result is somewhat counter-intuitive, as one would expect the driving force 
for absorption (and, hence, the absorption rate) to increase with increasing concentration of free chains. A slowing 
down of absorption kinetics with growing size L and concentration Ci of the free chains has been experimentally 
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observed 62j in a porous medium (activated carbon) which resembles in certain aspects the polymer brush. In Fig.[8jD 
we show the variation of the absorbed amount, T(t) oc , for the critical attraction e po — 1.70 - see HVBl We point 
out that this well pronounced power law increase of T was observed only at this particular value of e po whereas for 
e po = 2.00 where most of our kinetic measurement were performed, no simple T — t relationship was found - cf. Fig. 
[7^,. Thus, in a sense, the particular kinetics of absorption underlines the special role of the critical compatibility 
between brush and free chains. 




FIG. 9: (a) Changing concentration profile <j> (z) of free chains with time elapsed after a quench from \ = to \ — 2.0 
from DFT data. The time is given in logarithmic coordinates. Here the mean concentration c = 0.1875 and the time unit 
corresponds to 25000 MCS. The polymer brush is located at < Z < 12.5. (b) Variation of the flux of free chains into the 
polymer brush with time for two concentrations c = 0.0625, 0.25 i.e., Mo = 8, 32, and L — 64. 

In order to shed further light on the observed behavior, Fig. |H] shows DDFT results for the time-dependent density 
profile <fi (z,t) (iV =24, left panel) and flux j(z,t) (N a =8 and 32, right panel). In the left panel, one observes two 
"ridges" in (f> (z,t) at all times - a principal ridge, initially located in the bulk above the brush, moves gradually 
inside, while another (smaller) ridge is located near the opposite (bare) wall and gradually disappears still moving in 
the bulk. From the right panel, one can see that for higher concentration (N a =32) the flux prevails over the lower 
concentration one inside the brush and at shorter times (thereby explaining higher initial values of T seen in Fig. [8] 
for larger values of N a ), while the situation is reversed outside the brush at longer times. The latter behavior is 
presumably due to higher mobility at lower concentrations and explains the decrease of slope a with N a seen in the 
inset of Fig. [8] 



V. DISCUSSION 



In this work we studied a scarcely explored yet important aspect of oligomer and linear macromolecule absorption 
in a polymer brush - the case of (more or less) good compatibility between species in the bulk and grafted chains. 
Starting from oligomers (mono- and dimers) and going up to chain lengths L which exceed twice the length N of the 
grafted chains, we have determined the conformation of the absorbed species, the absorbed amount T, and absorption 
kinetics (the propagation rate into the polymer brush) at different concentration of the free chains for two cases of 
moderately to very dense polymer brushes. In addition, by combining Monte Carlo simulations with DFT and SCFT 
calculations, we have substantially broadened the range of lengths of the grafted chains to 32 < N < 256 in order to 
test more comprehensively our findings. 

The most salient, and - to some extent - unexpected features of linear chain absorption in a polymer brush that we 
find are: 

• the dramatic increase in adsorbed amount T(L) with growing chain length L, and 

• the significant slowdown of absorption kinetics with growing concentration (i.e., with the increase of the starting 
gradient in density) of the free chains 

Besides these static and dynamic properties of polymer absorption in brushes, we find that both the absorbed macro- 
molecules and the brush itself largely retain their structure and conformation, as seen in quantities like R g , R e and the 
monomer density profile 4> p {z), for different length L and concentration <j) of the free chains, and different strength e po 
of attraction to the grafted chains. In particular, the degree to which the brush profile 4> p (z) is affected by absorption 
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is found to be much less that anticipated in some earlier theoretical predictions (42|. Nontheless, even within these 
small changes we observe a slight contraction of 4> p {z) at small absorbed amounts T while <fi p (z) gradually attains its 
extension roughly to that corresponding to zero concentration of free chains with growing T. 

An interesting finding which still needs deeper understanding is the observed existence of a critical compatibility 
X c < (i.e., brush-oligomer attraction e£ D ). At \ c we find both in MC as well as in DFT/SCFT that the energy of all 
absorbed species has a value independent of their size L whereas their entropy experiences a minimum. The critical 
attraction Cp is manifested by the existence of unique distance from the grafting plane where all monomer density 
profiles of the free chains intersect. Moreover, at Cp the kinetics of free chain absorption into the brush follows a 
clear cut power law with exponent /3 oc L 1 / 3 . Undoubtedly, much more work is needed until all these fascinating new 
features are fully understood. 

Last not least, we emphasize the finding of three distinct regimes in the kinetics of free chain absorption as far as 
the size of the free chains L is concerned. In the first regime the characteristic time for absorption r grows rapidly 

— 1/2 

with oligomer length L as long as the oligomer size R g oc L v ss a g remains smaller than the separation between 
grafting sites. The second regime is marked by a slower increase of r with L and ends roughly at L N. The third 
regime of absorption kinetics holds for L > N (i.e., the penetrating free chain cannot accommodate within the brush) 
and is characterized by a nearly constant t as far as length L is concerned. Interestingly, this rich kinetic behavior 
has been experimentally observed in absorption in porous media [(32| . 
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